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Abstract 


A possible means of optical sensing, based on a porous hyperbolic material which is inhltrated by a 
fluid containing an analyte to be sensed, was investigated theoretically. The sensing mechanism relies on the 
observation that extraordinary plane waves propagate in the infiltrated hyperbolic material only in directions 
enclosed by a cone aligned with the optic axis of the infiltrated hyperbolic material. The angle this cone 
subtends to the plane perpendicular to the optic axis is 0c- The sensitivity of 0c to changes in refractive 
index of the inhitrating huid, namely was explored; also considered were the permittivity parameters and 
porosity of the hyperbolic material, as well as the shape and size of its pores. Sensitivity was gauged by the 
derivative dOc/dub- In parametric numerical studies, values of dOc/dm in excess of 500 degrees per refractive 
index unit were computed, depending upon the constitutive parameters of the porous hyperbolic material 
and inhitrating huid, and the nature of the porosity. In particular, it was observed that exceeding large 
values of dOc/dub could be attained as the negative-valued eigenvalue of the inhltrated hyperbolic material 
approached zero. 
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1 Introduction 


Through judicious design, engineered composite materials — sometimes referred to as metamaterials — 
can offer unique opportunities for controlling the propagation of light. Hyperbolic metamaterials represent 
a particularly interesting category of such engineered composite materials [1]. The simplest example of 
a hyperbolic material, at least from a mathematical perspective, is provided by a nondissipative uniaxial 
dielectric material whose permittivity dyadic is indefinite (i.e., it possesses both positive- and negative-valued 
eigenvalues). The hyperbolic-material concept may be extended to dissipative uniaxial dielectric materials 
whose permittivity dyadics possess real parts which are indefinite. Such materials are promising for a host of 
potential applications, including as platforms for negative refraction (without concomitant negative-phase- 
velocity propagation) [2, 3], subwavelength imaging [4, 5, 6], radiative thermal energy transfer [7, 8], as 
analogs of curved spacetime [9, 10], and as diffraction gratings which can direct light into a large number 
of refraction channels [11], for examples. Hyperbolic metamaterials may be fabricated as metal/dielectric 
composites, with the components arranged as layered sheets [12, 13] or as arrays of nanowires [14]. A simpler 
method of fabrication may be achieved via the homogenization of a random distribution of electrically-small 
ellipsoidal particles [15]. In this context it is notable that hyperbolic materials are also found in nature in 
certain spectral regimes; for example, graphite in the ultraviolet regime [16]. Greater scope for hyperbolic 
metamaterials may be offered by dielectric-magnetic materials whose permittivity and permeability dyadics 
are both indefinite (or whose real parts are both indefinite) [17]. A major challenge in developing hyperbolic 
metamaterials is to overcome dissipative losses. The use of metals such as silver at visible wavelengths and/or 
active component materials may help to minimize such losses [18, 19]. 

Our interest here lies in a novel optical-sensing application of hyperbolic metamaterials. We consider a 
porous hyperbolic material which is infiltrated by a fluid containing an analyte to be sensed. The presence of 
the analyte-containing fluid changes the optical properties of the hyperbolic material in manner that may be 
harnessed for optical sensing. Specifically, extraordinary plane waves propagate in the infiltrated hyperbolic 
material only in directions enclosed by a cone aligned with the optic axis of the infiltrated hyperbolic material. 
The angle this cone subtends to the plane perpendicular to the optic axis, namely 9c, may be acutely sensitive 
to the refractive index of the infiltrating fluid. Thus, by monitoring the magnitude of 9c, the concentration of 
analyte in the fluid may be gauged. In the following we explore the sensitivity of 9c to changes in refractive 
index of the infiltrating fluid; also considered are the porosity and permittivity parameters of the hyperbolic 
material, as well as the shape and size of its pores. The theoretical approach adopted is based on the extended 
version [20] of the well-established Maxwell Garnett homogenization formalism [21, 22], which requires that 
the hyperbolic material is highly porous and that the pores are much smaller than the wavelengths under 
consideration. 
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As regards notation, in the following vectors are underlined, with the" symbol denoting a unit vector. 
Thus, the unit vectors aligned with the Cartesian axes are {i, y,l.}- Double underlining indicates a 3x3 
dyadic (i.e., a second-rank Cartesian tensor), with / = xx + yy + H being the identity dyadic. The 
permittivity and permeability of free space are written as Cq and /Xq, respectively; kg = uiy/Cofio is the free- 
space wave number, with uj being the angular frequency. The real and imaginary parts of complex-valued 
quantities are delivered by the operators Re{ • } and Im{ • }, respectively; and i = 


2 Plane waves in hyperbolic materials 


Let us consider a homogeneous uniaxial dielectric material characterized by the Tellegen constitutive relations 

[ 22 ] 

Hit) = £ * Kit) 1 

}, ( 1 ) 

B{r) = HgH{r) j 

where the permittivity dyadic 

e = £■'■ (I — 11) + (2) 

The optic axis of the material is parallel to the unit vector z [23]. Largely we focus on nondissipative 
materials for which the permittivity parameters £■’■,£11 G K., but this restriction is relaxed in §4.2 wherein 
weakly dissipative materials (for which £■’■,£11 G C) are considered. 


Now suppose that plane waves with field phasors 

Mit) = Ejq exp {ik • r) 

> 

Kit) = iLo ' t) 


(3) 


propagate in the hyperbolic material characterized by Eqs. (1) and (2). In the case of uniform plane waves, 
propagating in a nondissipative material, the wave vector components are given as 


k = kxi + kyy + kzZ, {kx,ky,kz &M.). 

The source-free Maxwell curl postulates 

V X E(r) — = Q 

> 

V X h _{ l ) + = Q 

combined with the constitutive relations (1) provides the vector Helmholtz equation 

[(VxX) . (Vxl)-^oa;2|] •^(r) = 0, 


(4) 


(5) 


( 6 ) 
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which yields the dispersion relation 


{kl + ky + kl- [kl + ky) e-*- + /cfe'l - e-'- 


= 0 


(7) 


for plane waves specified by Eqs. (3)i. Two distinct plane-wave solutions emerge from Eq. (7): the ordinary 
plane wave whose wave-vector components satisfy [24] 


e = 


and the extraordinary plane wave whose wave-vector components satisfy [24] 




Cn 


( 8 ) 


(9) 


Note that the ordinary and extraordinary plane waves can be readily distinguished from each other by their 
polarization states. For example, if propagation is restricted to the xz plane, then the vector Helmholtz 
equation (6) yields the eigenvector solutions 


^0 = 


Eyy 

Exz 


e-*- kx 


for the ordinary plane wave 
for the extraordinary plane wave 


( 10 ) 


where the amplitudes Ey and Exz may be determined from the boundary conditions. 

Conventionally, in the realm of crystal optics e is positive definite (i.e., e"'" > 0 and > 0); therein, the 
ordinary dispersion relation (8) has a spherical representation in {kx,ky,kz) space and the extraordinary 
dispersion relation (9) has a spheroidal representation in {kx,ky,kz) space [23]. However, our interest here 
is in the case where e is indefinite. To be specific, suppose that e"*" > 0 and < 0. In this case, the 
dispersion relations for the ordinary and extraordinary plane waves are dramatically different. The ordinary 
dispersion relation (8) has a spherical representation in {kx,ky,kz) space but the extraordinary dispersion 
relation (9) has a hyperbolic representation in {kx,ky,kz) space. Accordingly, materials characterized by 
such an indefinite permittivity dyadic are known as hyperbolic materials [1]. 

An obvious distinction between the ordinary and extraordinary plane waves for hyperbolic materials 
is that ordinary plane waves can propagate in all directions whereas extraordinary plane waves can only 
propagate in a restricted range of directions. Specifically, extraordinary plane wave directions are bounded 
by the cone 

= {kl + kl) tan^ 0^, (11) 

which subtends the angle 

.-1 


= tan 


( 12 ) 
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to the kx-ky plane. The limiting values of the extraordinary wave-vector components are asymptotic to 
this cone. By way of illustration, the extraordinary wave-vector components kx and k^ in the ky = 0 plane 
are plotted in Fig. 1 for the case where the permittivity parameters £■*■ = 3.5eo and = —2eo. For wave- 
vector directions lying outside those prescribed by the cone (11), the extraordinary dispersion relation (9) 
yields only evanescent (i.e., nonpropagating) solutions. That is, extraordinary plane waves with wave-vector 
component k^ = kcosip, where k = only propagate in the directions given by 0 < '0 < (tt/S) — 9c- 

The transition between propagating solutions and evanescent solutions that occurs at 0 = {it 12) — 9c has 
the potential to be harnessed for optical sensing. 


3 Homogenization 


As a platform for optical sensing, let us consider a mixture of two materials: material ‘a’ is a porous 
hyperbolic material as characterized in Eqs. (1) and (2); and material ‘b’ is a fluid of refractive index rif,, 
containing a to-be-sensed analyte. The permittivity dyadics of the two component materials are written as 


la = +eaii 


(13) 


lf, = ^4 J 

with eb = All material ‘a’ pores are assumed to have the same shape and orientation. Each pore is 

taken to be spheroidal in shape, with the rotational symmetry axis of each pore being aligned with the optic 
axis of material ‘a’. Thus, the surface of each component pore relative to its centroid is traced out by the 
vector 


Ls = P^' t 


(14) 


wherein p > 0 is a measure of the linear dimensions of the pore; the shape-and-orientation dyadic 


g = U-^{l-zz) + uHz, > o) , (15) 

and f is the (unit) vector that prescribes the surface of the unit sphere The pores are randomly distributed 
in material ‘a’ with volume fraction /„. Hence, the volume fraction of material ‘b’ is 1 — fa- A schematic 
diagram of the two component materials is provided in Eig. 2. The particular nanostructure of the hyperbolic 
host material (at length scales smaller than the pores it contains) could be engineered as a one-dimensional 
[14], two-dimensional [12, 13] or three-dimensional [15] structure; or indeed, in principle, this material could 
even be a natural material [16]. However, the particular nanostructure of the hyperbolic host material (at 
length scales smaller than the pores it contains) is not directly pertinent to the proceeding analysis. 


Provided that the size parameter p is much smaller than the wavelengths under consideration, the mixture 
of materials ‘a’ and ‘b’ may be regarded as being effectively homogeneous. The permittivity dyadic of the 
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corresponding homogenized composite material (HCM), namely [22] 


= hcm ^hcm {L 


(16) 


can be estimated using the well-established Maxwell Garnett formalism. For the case where material ‘a’ is 
highly porous, i.e., fa ^ 0.3, the Maxwell Garnett estimate is given by [22] 

-1 


= hem 


+ (l- faR^ 's) , 

wherein the polarizability density dyadic 

(la-ih) • [l + R’ (ia-ib)] 
Two depolarization dyadics appear in Eqs. (17) and (18). First, 

Rs 

which corresponds to a spherical Lorentzian cavity, with 

D° = —I 

= s 3ef,= 

being the depolarization contribution arising in the limit p —> 0 [ 22 ], and 


9eo 


(3 + 2ikopnb)L 


being the depolarization contribution arising from nonzero p [25]. Second, 

D = D° + D+ 


(17) 


(18) 


(19) 


( 20 ) 


( 21 ) 


( 22 ) 


which corresponds to a pore whose shape and orientation are specified by V_ per Eqs. (14) and (15), with 

{I-zz) + D°’Hz (23) 

being the depolarization contribution arising in the limit p —^ 0 [ 22 ], and 

Q+= {I-zz) + zz (24) 

being the depolarization contribution arising from nonzero p [25]. Herein the depolarization dyadic compo¬ 
nents [26, 27] 


DO’-L = 


1 


2et {U^r 


tanh y/1 — 7 


1 -H 


1 


1 


HO’II = 


VI - 7 V 7-iy 7-1 
tanh“^ VI — 7 ' 


( 7 -l)(C/ll)^ 


1 - 


VI -7 


(25) 
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and [ 20 ] 


D+’-L = -- 




8eo (C/J 


7 


Akopvi^b 


D+’W = -- 


ey 


4eo {uwy 


1-7 

^koP\ki 1 


(3 + 7 ) - 


1-7 


\/l — 7 tanh ^ a /1 — ' 


1 - 7 V 1 - 7 


1 — ^- % -^/l — 7 tanh ^ a /1 — 7 


(26) 


with and the dimensionless parameter 7 = {U'^/Uy'^ Usually, the standard 

version of the Maxwell Garnett formalism is implemented wherein the approximation D_^ = and ^ 
is made; i.e., the nonzero size of the component pores is neglected. Indeed, the standard version of the 
Maxwell Garnett formalism is implemented in the numerical studies presented in the following §4.1. However, 
in §4.2, the nonzero size of the component pores is explicitly taken into account. This requires the extended 
version of the Maxwell Garnett formalism wherein the depolarization dyadics and are as given in 
Eqs. (19) and (22), respectively [20]. Notice that the standard version is identical to the extended version 
implemented in the limit p —> 0 . 


4 Numerical studies 

We investigate the prospects for optical sensing in a fluid-infiltrated hyperbolic material by means of numer¬ 
ical studies based on the homogenization approach described in §3. The envisaged mechanism for sensing is 
delivered by the transition between propagating solutions and evanescent solutions that emerges from corre¬ 
sponding dispersion relation for extraordinary plane waves. The transition occurs when the angle between 
the wave-vector direction and the optic axis of the infiltrated hyperbolic material is ( 7 r/ 2 ) — Oc- Thus, in the 
following attention is focussed on the critical angle 6^ and its dependency on the refractive index rib of the 
fluid containing the analyte to be sensed. The derivative dOc/dub provides a measure of the sensitivity for 
the envisaged method of optical sensing. 

Two qualitatively different scenarios are considered. In §4.1, the standard Maxwell Garnett formalism is 
used to investigate the effects of porosity, pore shape and permittivity parameters of the porous hyperbolic 
material, as well as the refractive index of the infiltrating fluid, on the sensitivity measure dOc/drib. In 
§4.2, the extended Maxwell Garnett formalism is used to investigate the effect of pore size on the sensitivity 
measure dOc/drib- The use of the extended Maxwell Garnett formalism results in a dissipative HGM, as may 
be appreciated from Eqs. (26). The dissipation introduced by such an extended homogenization formalism 
may be attributed to coherent scattering losses arising from the nonzero size of the pores. This phenomenon is 
a general feature of extended homogenization formalisms. Strictly, the corresponding HGM is not hyperbolic 
in the sense defined in §2 because its wave-vector components are complex-valued. However, provided 
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that I Im | <C | Re |, the real parts of the extraordinary wave-vector components satisfy the 

dispersion relation 

Re + Re {ky}'^ Re{fc^}^ _ 
and the HCM may be regarded as being approximately hyperbolic. 

4.1 Standard Maxwell Garnett formalism 

Suppose that the porous hyperbolic component material is specified by the permittivity parameters = 
3.5eo and eH = —2eo and the infiltrating fluid by the refractive index ni, = 1.3. The corresponding estimates 
of the HCM’s permittivity parameters, as yielded by the standard Maxwell Garnett formalism, are plotted 
against volume fraction fa in Fig. 3 with the relative shape parameter (C/H/G^) € {0.75,1.0,1.25}. In 
conformity with Eq. (17), the volume fraction range is restricted to 0 < /„ < 0.3. As fa increases, the HCM 
permittivity parameter increases approximately linearly from its starting point at 1.3^eo, and the rate of 
increase is largely independent of the value of /U^. In contrast, the HCM permittivity parameter is 
highly sensitive to the value of /U-^. As fa increases, makes the transition from being positive-valued 
to being negative-valued — which corresponds the HCM making the transition from being non-hyperbolic 
to being hyperbolic — at fa = 0.07,0.19 and 0.27 for = 0.75,1.0 and 1.25, respectively. 

The normalized wave-vector components of extraordinary plane waves are illustrated in Fig. 4 for the 
HCM whose permittivity parameters are displayed in Fig. 3. Only the volume fractions fa € (0.1,0.2,0.3} 
for = 1-0 are represented in Fig. 4. At fa = 0.3, the range of directions in which extraordinary 

plane waves can propagate is relatively large, with the critical angle being 0c = 62.7°. At fa = 0.2, this 
range of directions is much reduced with the critical angle being 9c = 78.8°. And at fa = 0.1 the material is 
not hyperbolic at all — the corresponding dispersion relation has a spheroidal representation in {k^, ky,kz) 
space. Parenthetically, the transition from spheroidal (or elliptical) to hyperbolic dispersion relations may 
be characterized by a sharp increase in photon density of states [28], which can result in increased rates of 
spontaneous emission from nearby emitters [29, 30]. Furthermore, for graphene-based hyperbolic materials, 
for example, this transition may be controlled by electrostatic biasing [31]. 

The derivative d9c/drn,, which is adopted here as a measure of sensitivity for optical sensing, is plotted 
against volume fraction fa in Fig. 5, for the same component materials as were considered in Fig. 3, with the 
relative shape parameter £ {0.75,1.0,1.25}. The presented ranges for fa are restricted to those 

ranges for which the HCM is hyperbolic; i.e., fa > 0.07 for (C/H/t/-*-) = 0.75, fa > 0.19 for (Ijll/t/-*-) = 1.0, 
and fa > 0.27 for {U^^/U-^) = 1.25. For all values of the sensitivity measure dOc/drn, increases 

rapidly as fa decreases towards its lower bound beyond which the HCM is no longer hyperbolic; and in this 




region of rapid increase, dOc/dri}, attains values in excess of 500 degrees per refractive index unit (RIU). 
Away from the region of rapid increase, for example at fa = 0.3 for (l/H /U-^) G {0.75,1.0}, dOc/duh attains 
values in the range 60-80 degrees per RIU. 

The effect of the pore shape on the sensitivity for optical sensing is explored further in Fig 6, wherein 
the derivative dBc/dnt is plotted against the relative shape parameter /[/■*' for the volume fraction fa € 
{0.2,0.25,0.3}. As in Fig. 3, the porous hyperbolic material is specihed by the permittivity parameters 
e}- = 3.5eo and eH = —2eo and the infiltrating fluid by the refractive index n/, = 1.3. We see that the 
sensitivity measure dOc/dni, increases uniformly, at a fixed value of porosity, as the pores become more 
elongated in the direction of the optic axis of the hyperbolic HCM. This general trend is independent of fa, 
but generally higher values of dOc/dni, are attained when fa is smaller. 

Next we turn to the refractive index Uf, of the analyte-containing fluid that infiltrates the porous hyper¬ 
bolic material. Its bearing upon the sensitivity measure dOc/drib is illustrated by Fig. 7, wherein dOc/drib is 
plotted against Ub for the volume fraction fa G {0.2,0.25,0.3}. The pores are assumed to be spherical for 
these plots; i.e., the relative shape parameter = 1.0. The permittivity parameters of component 

material ‘a’ are the same as for Fig. 3. In Fig. 7, dOa/dub increases uniformly, at a fixed value of porosity, 
as the refractive index Ub increases. Furthermore, the values of dO^/dub become exceedingly large as Ub 
increases towards its limiting value beyond which the HCM is no longer hyperbolic. As the volume fraction 
increases, this limiting value of Ub occurs at higher values of Ub- 

The influence of the permittivity parameters of the porous hyperbolic material upon the sensitivity 
measure dOc/dub is delineated in Fig. 8. Therein, dOc/dub is plotted against /cq and eH/eo for the volume 
fraction fa = 0.3, relative shape parameter (C/H/!/■*■) = 0.8, and refractive index of the analyte-containing 
fluid Ub = 1.3. The sensitivity measure decreases slightly as increases from Cq to 4eo, with the gradient of 
decent being greater at larger values of eL On the other hand, the sensitivity measure increases markedly 
as increases from —2eo to —1.3eo, with the gradient of ascent being greater at smaller values of ef. Thus, 
the greatest values of dOc/dub are attained at the smallest values of ef; largest values of 

4.2 Extended Maxwell Garnett formalism 

Now we consider the influence of pore size upon the sensitivity measure dOc/dub- For this purpose, the 
extended Maxwell Garnett formalism is required. The real and imaginary parts of the HCM’s permittivity 
parameters, as estimated by the extended Maxwell Garnett formalism, are plotted against the normalized 
size parameter kop in Fig. 9, with the relative shape parameter = 1.0 and volume fraction fa G 

{0.2,0.25,0.3}. As in Fig. 3, the porous hyperbolic material is specihed by the permittivity parameters 
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= 3.5eo and eH = —2eo and the infiltrating fluid by the refractive index Uh = 1.3. The real part of 
increases slightly as the pore size increases, with the values of Re being only slightly influenced by 

the volume fraction. In contrast, the real part of increases markedly as the pore size increases, with the 
values of Re being considerably influenced by the volume fraction. The imaginary parts of 

ej^cm both increase from zero as the size parameter increases, with the rate of increase being substantially 
greater for Im than for Im ■ These imaginary parts are only marginally effected by the volume 

fraction. 


In order to implement the sensitivity measure dOc/drih for the p > 0 regime, it is important to take into 
account the complex-valued nature of the HCM’s permittivity parameters. In regions where Re < 0 

and Re > 0, the magnitudes of the imaginary parts of and are both generally very small 

in comparison to the magnitudes of the real parts of e|cm) presented in Fig. 9. Exceptionally, 

when Re extremely close to zero (and is negative valued), Im can become relatively large 

compared to Re |e|cm} consequently the characterization of the HCM’s extraordinary dispersion rela¬ 
tion per Eq. (27) is not appropriate. The regime where Ri 0 is discussed further in Section 5. 

Therefore, excluding the regime where ~ 0, the HCM may be regarded as being approximately 

hyperbolic with the real parts of the extraordinary wave-vector components satisfying the dispersion rela¬ 
tion (27). Accordingly, here we extend the definition (12) of the critical angle, which marks the transition 
from propagating solutions to evanescent solutions emerging from the corresponding dispersion relation for 
extraordinary plane waves, to 


9c = tan ^ 






(28) 


That is, extraordinary plane waves propagate in the approximately hyperbolic HCM for directions given by 
0 < tp < (7r/2) — 9c, where ip = cos“^ (Re {kz} /Re {fc}) and the critical angle 9c is as given in Eq. (28). 


The derivative d9c/drii, is plotted against the normalized size parameter k^p in Eig. 10. Therein the volume 
fraction fa € {0.2,0.25,0.3}, the relative shape parameter (I/H/t/-*-) = 1.0 and the constitutive parameters 
of the component materials ‘a’ and ‘b’ are the same as for Eig. 9. For all volume fractions considered, the 
sensitivity increases uniformly as the pore size increases. The rate of increase of d9c/dni, is relatively modest 
at low values of p but increases rapidly as p increases. At a fixed value of p, the greatest values of d9c/drib 
are attained at the smallest values of fa- 
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5 Discussion 


A novel means of optical sensing has been explored theoretically. It is based on a porous hyperbolic material 
which is infiltrated by a fluid containing an analyte to be sensed. Extraordinary plane waves can only 
propagate in the infiltrated hyperbolic material in directions enclosed by a cone which subtends the critical 
angle O,. to the plane perpendicular to the optic axis of the hyperbolic material. Numerical studies utilizing 
the Maxwell Garnett homogenization formalism have revealed that the critical angle 9^ is highly sensitive 
to the refractive index of the analyte-containing fluid, and this sensitivity can be strongly influenced by the 
permittivity parameters and porosity of the hyperbolic material, as well as the shape and size of its pores. 
Specifically, the sensitivity increases as: (i) the refractive index of the fluid increases; (ii) the negative¬ 
valued eigenvalue of the hyperbolic material’s permittivity dyadic increases; (iii) the porosity increases; (iv) 
the pores become more elongated in the direction of optic axis; and (v) the pores increase in size. 

The preceding analysis is based on the well-established Maxwell Garnett homogenization formalism 
[21, 22]. The key assumptions underlying this approach are that the porosity of the hyperbolic host material 
is relatively high (i.e., fa ~ 0.3) and that the pores are much smaller than the wavelengths involved. Like 
most conventional approaches to homogenization, this approach is limited to materials which are spatially 
local. If the hyperbolic host material were to possess a characteristic length scale which was similar to 
the wavelengths involved, then the effects of spatial nonlocality may become significant and an alternative 
homogenization approach may be needed [32]. However, if the hyperbolic host material were envisaged to 
be a homogenized composite material arising from components made up of particles much smaller than the 
wavelengths involved [15], for example, then the effects of spatial nonlocality are likely to be negligible. 

It is notable that exceeding high sensitivities, as gauged by the derivative dOc/drib, can be achieved. For 
a HGM with permittivity parameter > 0 investigated in §4), the largest values of dOc/dub are 

attained when the negative-valued HGM permittivity parameter is very close to zero. Indeed, in the 
limit —>■ 0“, the critical angle 9c 7r/2 and the derivative d9c/dnb becomes unbounded. However, 

when losses are taken into account, in the immediate vicinity of Re ~ magnitude of Re 

can become relatively small compared to that of Im consequently the characterization of the HCM’s 

extraordinary dispersion relation per Eq. (27) may be inappropriate. But even for values of e^lcm away 
from the immediate vicinity of Re |e|cm} ~ magnitudes of the sensitivity measure d9c/dnb presented 

in Figs. 5, 6, 7, 8 and 10 compare favourably to sensitivities associated with optical sensing approaches 
based on surface-plasmon-polariton waves [33, 34, 35, 36, 37], Voigt waves [38, 39], as well as generic porous 
dielectric materials [40]. We note that the issue of relatively large losses that arises in the parameter regime 
where Re | | « 0 could, in principle, be addressed by the use of active component materials [41]. 
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The numerical results reported here bode well for the prospects of harnessing porous hyperbolic materials 
as platforms for optical sensing, and thereby adds to the already long list of potential applications for these 
exotic materials [1]. 
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Figure 1: Normalized wave-vector components kx/ko and k^/kg of extraordinary plane waves for the per¬ 
mittivity parameters e'*" = 3.5eo and = — 2eo. The ky = 0 plane is represented, for which the wave-vector 
components are bounded by the cone k^ = k^ tan^ 0c, where 9c = tan“^ 
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Figure 2: Schematic diagram showing the composite material to be homogenized using the Maxwell Garnett 
formalism. A highly-porous hyperbolic host material (labelled ‘a’) contains spheroidal pores which are filled 
with a fluid (labelled ‘b’). The electromagnetically-small pores all have the same orientation but they are 
randomly distributed. 
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Figure 3: Maxwell Garnett estimates of e^cm/^o and ej^^^/eo plotted versus volume fraction fa G (0,0.3) for 
the relative shape parameter (t/H/C/-*-) = 0.75 (blue, dashed curves), 1.0 (red, solid curves) and 1.25 (green, 
broken dashed curves), refractive index rif, = 1.3, and permittivity parameters e-f = 3.5eo and eH = —2eo. 
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Figure 4: Normalized wave-vector components of extraordinary plane waves for volume fractions fa £ 
{0.1,0.2,0.3}, with the relative shape parameter /U'^) = 1, refractive index Ub = 1.3, and permittivity 
parameters = 3.5eo and eH = — Scq. 
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Figure 5: The derivative dOc/dni, (in degrees per RIU) plotted versus volume fraction fa for the relative 
shape parameter (C/H/C/-*-) = 0.75 (blue, dashed curve), 1.0 (red, solid curve) and 1.25 (green, broken dashed 
curve), refractive index nt = 1.3, and permittivity parameters = 3.5eo and eH = —2eo. 
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Figure 6: The derivative dOc/dub (in degrees per RIU) plotted versus relative shape parameter for 

the volume fraction fa = 0.3 (red, solid curve), 0.25 (blue, dashed curve), and 0.2 (green, broken dashed 
curve), refractive index nt = 1.3, and permittivity parameters = 3.5eo and eH = —2eo. 
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Figure 7: The derivative dOc/dni, (in degrees per RIU) plotted versus refractive index nt for the volume 
fraction fa = 0.3 (red, solid curve), 0.25 (blue, dashed curve), and 0.2 (green, broken dashed curve), relative 
shape parameter ([/ll/C/“*~) = 1, and permittivity parameters = 3.5eo and eH = —2eo. 



Figure 8: The derivative dOc/duh (in degrees per RIU) plotted versus relative permittivity parameters 
f-a/i-o and eH/eo, for the volume fraction fa = 0.3, refractive index Uf, = 1.3, and relative shape param¬ 
eter (f/ll/t/-L) = 0.8. 
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Figure 9: Extended Maxwell Garnett estimates of the real and imaginary parts of ^Lm/^o 

plotted versus relative size parameter koP for the volume fraction fa = 0.3 (red, solid curve), 0.25 (blue, 
dashed curve), and 0.2 (green, broken dashed curve), relative shape parameter (C/H/f/-*-) = 1.0, refractive 
index Uh = 1.3, and permittivity parameters = 3.5eo and eH = —2eo. 
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Figure 10: The derivative dOc/dni, (in degrees per RIU) plotted versus relative size parameter koP for the 
volume fraction /„ = 0.3 (red, solid curve), 0.25 (blue, dashed curve), and 0.2 (green, broken dashed curve), 
relative shape parameter = 1.0, refractive index nt = 1.3, and permittivity parameters = 3.5eo 

and Ca = — 2eo. 
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